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Abstract We study the Efimov effect in a harmonic oscillator in the hyperspherical formulation, and 
show how a reduced model allows for a description that is a generalization of the Efimov effect in free 
space and leads to results that are easily interpreted. Efimov physics may be observed by varying the 
value of the scattering length, since in the regime where the trimers have a mixed harmonic oscillator 
and Efimov character, the inelastic properties of these states are still manageable. The model also allows 
for the study of non-universal Efimov trimers by including the effective range scattering parameter. 
While we find that in a certain regime the effective range parameter can take over the role of the 
three-body parameter, interestingly, we obtain a numerical relationship between these two parameters 
different from what was found in other models. 
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1 Introduction 

In 1970, Efimov described a phenomenon occurring in a three-particle system of identical bosons, in 
which there is an accumulation of bound states near the point of zero energy and diverging scattering 
length of the interparticle interaction [1 . Efimov was inspired by the physics of the nucleus of the 
tritium atom. However, the experimental signature for the existence of such Efimov states was not 
observed in nuclear physics, but in a gas of ultracold cesium atoms [2] by studying three-body recom- 
bination. The reason that properties predicted for nuclear physics few-body systems can be derived 
from ultracold atoms lies in the fact that the two systems are connected through universal relations, 
that are irrespective of the specific nature of the interparticle interactions. 

Since the first observations in cesium, the Efimov effect has received a lot of attention [31 13] . Several 
groups now have observed Efimov physics by studying three-body recombination in a variety of atomic 
systems [SlEllT, 8; 9; 10, .11; 12; 13, 14J. Evidence for the existence of trimers appears as some distinct 
features (resonances) in the recombination rate coefficient. In Fig. [T] the energies of the Efimov trimers 
are shown, as a function of the inverse two-body scattering length 1/a. 

A crucial advantage of ultracold atoms with respect to nuclear systems is that the interparticle 
interactions can be easily manipulated using Feshbach resonances [151 116j , by varying the magnetic 
field. The two-body scattering length, which is obtained from the scattering phase shift d{k) as a = 
— limj._j.o tan(5(fc)/fc, diverges on resonance to infinity, and is typically described by the dispersive 
formula 
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Fig. 1 A schematic description of the Efimov effect: shown is the wavenumber K ~ VE, indicating dimer and 
Efimov trimer states as a function of the inverse scattering length 1/a. The horizontal line K = Q indicates 
the threshold to the three-body scattering continuum (gray). Simultaneously, the blue region marks the atom- 
dimer continuum for a > 0, and the green zone indicates the trimer region. The dimer binding energy scales 
as —1/a. Only a few Efimov states are shown. The figure is not to scale: for clarity reasons we plotted instead 

oi K the quantity //^^^sin^, and instead of 1/a the quantity _ff^''''cos^ (cf. Eq. ([3|l). 



with ahg the background value of the scattering length, B the magnetic field, Bq the resonance position, 
and AB the width of the resonance. This possibility to change the scattering length allowed for the 
observation of distinct Efimov features in the three-body recombination rate, connected to two different 
Efimov states in a row [6||8], confirming the scaling law [UllI] given by 

E'^+^/E^ = e-^'^/"" « 1/515, (2) 

with So — 1.00624 a constant. Moreover, recent developments extended the universal few-body physics 
to the domain of four-body states which was theoretically predicted [18- and experimentally verified [SJ 

US]. 

In most of the experiments, Efimov physics has been observed only indirectly via three-body re- 
combination [17 . Due to the short lifetime of these states, it is difficult to probe the Efimov states 
directly. These weakly-bound trimers are not only unstable in collisions with other particles, but they 
are also intrinsically unstable by itself, as trimers may decay to a deeper lying dimer state and an 
unbound atom. However, the group of Jochim in Heidelberg recently succeeded in associating Efimov 
trimers directly via a new and promising approach using radio- frequency fields |20j . 

Optical lattices might provide more experimental control, and potentially give rise to an enhanced 
lifetime of Efimov trimers. Efimov physics can be studied in this way by loading three atoms in each 
lattice site. When the lattice is deep enough, single sites act as harmonic oscillator potentials where 
individual trimers can be trapped. A possible advantage is then that the interparticle separation of 
the trimers can be increased to the scale of the harmonic oscillator length, while their properties are 
still determined by Efimov physics. In this way the overlap of these states with deeply bound dimers is 
reduced, thereby suppressing the process of three-body recombination. The influence of other limiting 
processes, where four or more particles involved, will be reduced as well. 

Jonsell, Heiselberg and Pethick were the first to introduce an harmonic oscillator term in a three- 
body system [3T], by adapting the approach of Ref. [25]. For diverging scattering length, they found 
exact expressions for the bound states of three identical bosons, while they obtained numerical results 
for finite scattering length. Later, Werner and Castin [53] obtained additional analytical results for 
diverging scattering length, and complimentary first-order corrections for large scattering length. Stoll 
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and Kohler [24] performed numerical simulations with a tunable separable potential to find the bound 
states of three particles in a trap. Th0gersen et al. found numerical solutions of three particles in a 
harmonic trap |25j , mainly to study the non- universal corrections to the Efimov effect that come from 
the effective range. More recently, Liu, Hu and Drummond derived an analytical description of three 
attractively interacting fermions in a harmonic trap [2^. 

A convenient description of the Efimov effect in free space is presented in Ref. jl7j . in which polar 
coordinates H and ^ are introduced, which are related to the wave number or energy parameter K 
and the scattering length a via 

K = Hsm.i, (3a) 

-=HcosS,. (3b) 
a 

Using Efimov's radial law, it is derived that 

for a function A : [— tt, — 7r/4] — > M. Here is the three-body parameter that defines the trimer 
energies exactly on resonance as EJj^ — — e^^"'^/*''/).^K^/m, with m the mass of a single atom, and is 
a result of the short-range behavior of three particles being close together. The combinations of K 
and ^ that yield bound states describe a family of curves that can be mapped on top of each other by 
multiplying the distance of the curves to the origin by a factor e^/'*", see Fig.[l] 

In this paper, we extend the description of Ref. |17| and show how by introducing the harmonic 
oscillator length as a new parameter, and by using spherical coordinates in parameter space, the scaling 
behavior in free space very naturally generalizes to a situation with three particles in a harmonic trap. 
In particular, the curves that are shown in Fig. [l] will be replaced by surfaces that are mapped on top 
of each other. In the derivation, we use an approach slightly different from earlier work, that calculates 
the equivalent of the function A by solving an ordinary differential equation. Additionally, we also 
include the coupling terms, we investigate the stability of the trimers in the lattice, and we study the 
non-universal corrections coming from the effective range parameter. 



2 Framework 

We consider a system of three bosons with pair interactions trapped in a harmonic potential. The 
corresponding Hamiltonian is 



H 



where Oho = \/h/mu} is the harmonic oscillator length, uj the trap frequency, and Vi is the interaction 
between the particles j, fc, with (i, j, fc) a cyclic permutation of (1, 2, 3). For convenience, the Hamilto- 
nian is scaled by a factor m/fi?. Throughout this paper, we will denote the eigenvalues of H by E, and 
the corresponding trimer energies by Et = h^E/m. To find the bound state energies, we closely follow 
the method described in Ref. [22] and Ref. [T7] . We denote the permutation operator corresponding to 
one cyclic permutation of (1, 2,3) by P and transform to the Faddeev operator 



F 



\Y,A, + -^Y.^j + V{I + P + P^). (6) 

^"ho 



If ^ is an eigenfunction of F, (/ + _P + P^)W is an eigenfunction of H , provided it is nonzero. In case 
{I + P + P"^)^ = Q,^ IS called a spurious solution. We separate the center of mass motion, and consider 
the relative system in hyperspherical coordinates. For that, we introduce Jacobi-coordinates according 
to 

r^Jfc = y|(r,- i(rj +rfe)), (7) 
T^3k = J\{rj-Vk). (8) 
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The hyperradius R and the hyperangle a are defined by 



jfe; 



i?COS Ui, 

Rsin Ui. 



(9) 
(10) 
(11) 



The remaining hyperspherical coordinates are given by the unit vectors ii^jk and ijk- 

We restrict ourselves to states with zero total angular momentum. The advantage of considering 
the Faddeev rather than the Schrodinger operator is that subsystem angular momentum comes in 
only at second order in the short-range potential [27l. This motivates us to neglect the subsystem 
angular momentum, by averaging over the Jacobi coordinates f ijfc and Vjk- In this approximation, the 
Schrodinger wavefunction tf'(ri, r2, rs) is given in terms of reduced Faddeev wavefunctions il^iR, a): 

if(ri,r2,r3) =V^(i?,ai) + 7/;(i?,a2)+V^(i?,a3). (12) 
In the differential Faddeev equation [ 17] we then make an adiabatic hyperspherical expansion 

1 



^(i?, a) — 



27ri?5/2 sin(2a; 



Y,UR)K{R.a). 



(13) 



For every value i?, which is treated as a parameter, the functions 0„(i?, a) of hyperangle a form a 
complete set, and are solutions to the eigenvalue equation on L^((0,7r/2)) 



Xn{R) 



sin a) 



i{R,a) 



^Tr-\^■!T-a\ 



3 J\^^T-a\ 



4>n{.R, o/)da' 



, (14) 



with Dirichlet boundary conditions 4>{R,0) = (j){R, n/2) = 0. The eigenvalues A„(i?) depend paramet- 
rically on R. The set of hyperangular functions (j}n{R,Oi) are generally not orthogonal functions in a. 
Therefore, we define 



Gnm{R) 



0* (i?, a)(t>rn{R, a)da, 



(15) 



which allows us to derive the following coupled set of differential equations for the functions /„ 



4i?2 J 



1 1 

2<: 



R' 



A„(fl)-4 
2i?2 



.fn{R)- 



E 



2P 

£j1 ri. 



^{R)Q^+Qam{R) 



U{R)^EUR). (16) 



The coupling terms in these equations are defined by 



k 



i'k{R^OL)-Q^(l^m{R-,a)da, 



K{R^0i)-Q^4>ni{R: a)da. 



(17a) 
(17b) 



This framework is very similar to the work of Refs. |221 117) . and the harmonic oscillator term 
is taken into account as in Ref. |21j . For now, we neglect the coupling terms conform the adiabatic 
hyperspherical approximation. However, later on we wil study the coupling between different channels 
more carefully, as the harmonic oscillator gives rise to degeneracies between the Efimov-like solutions 
in the zeroth channel and the continuum-like solutions of the first and higher channels. 
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Hyperangular problem 

The effect of the interaction potential V is now encoded in the eigenf unctions A„(i?). Parametric 
solutions for A„ can be found by fixing the parameter R, while solving for the angular problem. The 



functions A„(i?) will then serve as input for the radial equation (16), and when the coupling terms are 



ignored, solutions are obtained by solving an ordinary differential equation (ODE). 



To find an approximation of A„, one can study solutions ij^iR, a) of Eq. ( 14) for small values of a 
and values of a close enough to 7r/2. Matching of the solutions at intermediate values of a then yields 
an approximation of A„. This procedure has been described in Ref. [22^. The equation found in this 
manner is 



cos 



[Xl/^) ^A-i/^sin (Ay^l) = -V2A-i/2sin (^Xl/^) R kncomkn)), (18) 



where kn = ^JXnjlR? . Note that A„ = 16 always solves this equation, but it gives rise to a spurious 



solution. The right hand side of equation ( 18 ) can be conveniently expanded using the effective range 



approximation for the s-wave scattering phase shift of the two-body problem 

/cflC0t^(fc;^)--- + ir,4 + 0(4), (19) 
a I 

with Tg the effective range parameter. When this parameter is neglected, we find the approximating 
equation for A„ = A„(i?/a) that is at the basis of the derivation in Ref. [22}: 



cos 



- ^A-V2 3in (aV^i) = ^A-V2(^),i, (aV^I ) ^. (20) 



By taking into account the effective range as well, one obtains an approximation that is more accurate 
with respect to the two-body physics near narrow Feshbach resonances. This has been proposed in 
Ref. [21] as a regularization procedure and has been readdressed in Ref. [53]. We will come back to 
this issue in Section |4j 

Hyperradial problem 

With the approximation of the channel eigenvalues at hand, we may turn to the hyperradial prob- 
lem. In the adiabatic hyperspherical approximation, we ignore the coupling between the different 
hyperradial channels and only take the scattering length into account in the approximation of A„. In 
this setting, the Efimov effect was derived in Ref. [H], and repeating it in this context allows for a 
comparison in the limit of Cho ^ oo- The uncoupled system of equations is given by 



XnjR) - I R^ 

' dR^ ^ i?2 + 

on {Rq,oo), with a boundary condition 



fn{R) = 2Efn{R), (21) 



/„(i?)^0, i?^oo, (22) 

and an appropriate boundary condition on the edge of the inner region Rq where the short-range 
potentials become important, which is given in terms of the three-body parameter k^. The exact 
boundary condition will be specified later. 

We first concentrate on the zeroth channel, i.e. n = 0, which is the only channel that has an 
attractive potential, and therefore the only channel in which the Efimov effect is found. We introduce 
spherical coordinates (T, d, ^) according to 

/s: Tsini^sinf, (23a) 

-=rsini?cosC, (23b) 
a 

— =rcosi?. (23c) 

Oho 
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and write down the differential equation for /(cr) — fo^a/T) as 

/(a) =2sgn(Osin2i?sin2e/(a), 



9^ Aq (cr sin 1? cos ^) 



cos ■& 



witfi long-range boundary condition 



0, cr 



(24) 



(25) 



For sma ll values of cr, the model does not provide us with a natural boundary condition. Solutions of 
Eq. (24) have a log-periodic behavior as cr —)■ 0, given by /(cr) ~ Y^sin(so log(Kcr/r) -|- ao) for some 
value of K. Here sq = Ao(0), and 



ao = --So log 2- T^Arg— — . 

2 2 1(1 — iso) 



(26) 



Fortunately, the real short-range potentials Vi provide us with a natural boundary condition at Rq 
encoded hy n — k^,, where is the three-body parameter that accounts for the short-range physics, 
which was mentioned already in the introduction. This gives rise to a short-range boundary condition 
at Ro'. 



-r, — r ^ n ^ cot(so log(K*i?o) -I- ao). 

/(cro) 2cro ao 



(27) 




Fig. 2 Surfaces that represent trimer bound states in the lowest channel, as a function of a ^, K and a^^^ 
(blue) . Only three surfaces are shown, but in fact there is an infinite number of them, which all can be mapped 
on top of each other by multiplying the distance to the origin by powers of 6^^"°. These surfaces are shced by 
a plane (red transparant), for a given aho that corresponds to the height of the plane. The intersection lines 
(red) represent bound state solutions for combinations of the scattering length a and the wave vector K, for 
this value of Sho- The picture is not too scale: we used similar scaling relations as used in Fig. [l]in order to 
show the surfaces more clearly. 



We may determine for which parameter combinations a solution exists by a shooting method from 
large values of cr. A particular convenient method makes use of a certain phase variable 4>. One can set 
up an equation for that is much alike the variable phase equation, but is adapted to fit the boundary 
condition at short distances very well. Henceforth, it will play a role similar to the function A as we 
mentioned in the introduction. We introduce the transformation 



/'(a) = {g{a) + h{a) cot(0(a) + f{a). 



(28) 
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Fig. 3 The trimer bound state energies of cesium in the zeroth and first hyperspherical channels. The result is 
obtained for a trap frequency lj = 1000 s~^. A deeper bound state, with much lower energy, can not be shown 
in this graph. Also indicated is the dimer binding energy (solid blue line) 



The functions g, h > and ip may be chosen to our convenience. Continuity properties in the parameters 
are easier to derive considering (j) than when one considers / directly. A particular choice 

/i=^, (30) 

(T 

^/j = sologcr + ao, (31) 



is highly compatible with the boundary condition (27 1. The equation for becomes 
'Ao(crsin-(?cos^) + Sq 



2a sgn(0 sin^ i9 sin^ ^ + cos"* da^^ sin^ (sq log a + uq + (f)ia)) = 0, (32) 



which is a nonlinear first order differential equation. Given values for ^ and we denote by (f>{(T; ^, i?) 
the solution to (32 1 that satisfies ^(ct) + sologcr + ao t 0. With this definition, the problem ([2| has a 
solution if and only if 

0(Ti?o; ^, I?) - so log(K*/r) mod tt. (33) 
Consequently, eigensolutions only exist if there exists an n G Z such that 

T = K,e-'^(™°'«''')/"''e-"''/"«. (34) 

We would like to emphasize here the similarity with the description of the Efimov effect in free space 
given by equation Using the asymptotic behavior of Aq for small arguments, one can easily show 
that (/'(cr) converges to a finite value for a 10. The function 0(0; ^,z?) can be shown to be continuous 
over the region. For -0 = and ^ G [— tt, — 7r/4], the Efimov effect in free space is retained. Hence, 
0(0;^, 0) is an approximation of — Z\(f)/2. 

Trimer solutions in a harmonic potential 

The trimer solutions are most easily represented by surfaces in three-dimensional (a"^ , K,a^^) 
space, as defined in Eq. (23). In Fig. [2] such surfaces are shown for the lowest channel, which can be 
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Fig. 4 Crossings in Fig. |3| are replaced by avoided crossings when the coupling between the zeroth and first 
hyperspherical channels is taken into account. 



mapped on top of each other by multiplication of the distance to the origin by powers of e^l^^ . Fig. [2] 
can be regarded as a generalization of the two-dimensional Fig. [T] When these surfaces are sliced with 
a plane corresponding to a fixed harmonic oscillator length aho, we obtain curves that express wave 
vectors K in terms of the scattering length a, and an example can be seen in Fig. l3| When taking 
the limit of ^ 00, the Efimov effect in free space is obtained as depicted in Fig.^ In the region 
^ e [— TT, — 7r/4], the intersections of the surfaces will tend to the curves that describe the Efimov states. 
In the region ^ g (— 7r/4, tt), the energy levels come closer and closer together, expressing the fact that 
in the limit, this region is filled with continuum states. The intersection of the surfaces with the a^^ = 
plane precisely corresponds to the exact results obtained by Jonsell et al. and Werner [5TJ [30] • 

As an example, we show in Fig. [3] the trimer binding energies of cesium as a function of scattering 
length with a fixed oscillator length, for the lowest and next to lowest hyperspherical channels. The 
parameters are k* = 1.93 x 10~'*aQ ^ [2], with gq the Bohr radius, Rq — lOOao, and Oho = ISOOOao, 
which corresponds to a trap frequency of approximately 1000 s~^. 

It is interesting to note the change of energy scales in this figure. Below the thresholds of the 
trimer and atom-dimer continuum (cf. Fig. [T]), which would exist for the free-space Efimov effect, the 
bound states show the characteristic scaling law of Eq. [2] Above these thresholds, the bound states 
show a spacing typical for the harmonic oscillator potential which is approximately 2fi,u;. This change 
of exponential scaling to linear scaling of the trimer energies makes that the accumulation of bound 
states in the point of zero energy and infinite scattering length disappears. Effectively, the Efimov 
trimers are "pushed" out of the trimer region through the thresholds and aquire a strong harmonic 
oscillator state character. 

There is also a clear difference between the zeroth and first hyperspherical channels. In this respect, 
it is interesting to note the different limiting behavior of the solutions for small positive and negative 
scattering length. For small and negative scattering length, there is no weakly-bound dimer state 
in the interatomic potential, and the solutions have the character of three weakly-interacting atoms 
in a harmonic oscillator potential. For small and positive scattering length, there is a weakly-bound 
dimer state, and for the zeroth hyperspherical channel the solutions have the character of a weakly- 
interacting dimer and an atom in a harmonic oscillator potential. These regimes are connected through 
solutions which have a strong Efimov character. However, for the first hyperspherical channel the 
situation is different. Here the solutions behave more as those of two particles undergoing a Feshbach 
resonance |31i l32j , which as a result interpolate between two different harmonic oscillator levels. 
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Fig. 5 The width Ft of cesium Efimov states in a trap, with u = 1000 s~^ and rj, = 0.06, for the first (red), 
second (black), third (green) and fourth (blue) trimer states depicted in Fig.js] Also, the width of a free Efimov 
trimer corresponding to the first state is shown (black dashed). 



It should be possible to investigate this change of character experimentally, by using similar tech- 
niques as have been done by Stoferle et al. [32]. Here the creation of dimers was investigated in an 
optical lattice while sweeping the magnetic field through a Feshbach resonance. For the trimers, such 
adiabatic sweeps through trimer solutions should allow for a creation of states with a character changing 
between harmonic oscillator-, Efimov- and dimer-|-atom- type solutions. 

So far, we have applied the adiabatic hyperspherical approximation, in which we completely ignored 
the coupling between the various hyperspherical channels. While this is generally a good approximation, 
we see from Fig. [3] that there are degeneracies between the zeroth and first channel solutions, which 
should be lifted when the coupling is taking into account. In order to demonstrate this, we introduce 
an approximation of the coupling. We calculate the coupling terms Pmn and Qmn of Eq. (17 1 using 
first order approximations of the functions and find in this way combined solutions which involve 
both channels. In free-space, the coupling between these two hyperspherical channels is directly related 
to the process of three-body recombination, where three colliding particles relax into a weakly-bound 
dimer and a third particle, with a surplus of kinetic energy |17| . Interestingly, for our lattice solutions 
we observe that this coupling transforms the crossings in Fig. [3] into avoided crossings, as shown in 
Fig. [4j This could be experimentally interesting, since using magnetic field sweeps would allow for 
adiabatic transitions from one hyperspherical channel to an other. 



3 Stability of trimers in a harmonic potential 

The change of character of the trimer solutions, which was discussed above, will have a profound 
impact on the stability of these states. So far, we implicitly assumed that the trimers in the harmonic 
oscillator potential are stable. However, when the three particles are close together the trimer can 
dissociate into a deeply-bound dimer and a third atom, with large kinetic energy. Here, we investigate 
the stability of the trimers by applying the method of Braaten and Hammer |17j . that regards the 
trimers as scattering resonances in a dimer-atom scattering process. In this respect, the trimer energy 
Et acquires a width Ft, where the latter can be used to calculate the lifetime h/Fr of the trimers. In 
this method, an inelasticity parameter ry* is introduced to account for the decay into various deeply 
bound states, which can be regarded as the imaginary part of a complex phase factor that describes 
the reflection of hyperradial waves in the inner region. Effectively, this comes down to changing the 
boundary condition at to 

= TTB- + -5- cot (solog(K*i?o) + ao + "/*) • (35) 



10 



For given values for the scattering length and the harmonic oscillator length, we are then able to find 
the associated complex eigenvalues 

where the trimer energy Et < and the width Ft > 0. 

As an example, we calculated the linewidths or cesium, with a trapping frequency to = 1000 s""'^, 
K* — 1.93 X 10~^ao, and rj^, = 0.06 [5]. The short-range boundary condition was implemented at a 
value Rq = lOOflQ. For some of the Efimov trimers shown in Fig. [3j the linewidths are shown in Fig. [5j 
The change of character of the trimers can be clearly seen from this graph: for values a/aho ^ —2 the 
states have a harmonic oscillator character, and the widths are relatively small. However, the Efimov 
effect shows up via a strong increase of the widths towards a/oho ^ —2, which is the parameter region 
where in free-space a new Efimov state would enter the tree-body system. It is also the region where 
a local peak in the three-body recombination rate would occur (trimer resonance). Here, we see that 
at this point the width of the lowest trimer keeps increasing, while for the higher trimer states the 
width saturates and displays an oscilating behavior. This is not surprising, since the latter trimers are 
more similar to three-particle continuum states. For comparison, the linewidth of the free-space Efimov 
trimer compared to the lowest lattice trimer state is shown in the region where it exists. These lines 
overlap mostly, and can be more than two orders of magnitude above the widths for the higher trimer 
states. At the point where the trimer energies approach the atom-dimer threshold, around a/aho = 34, 
the widths go down again to the level of the higher trimer states. The width of the lattice trimer is 
there somewhat larger than the free-space trimer, probably due to the additional confinement of the 
weakly-bound dimer state caused by the lattice. 

While the Efimov effect shows up in the excited states by a clear deviation from the harmonic 
oscillator levels, these states benefit from a much smaller linewidth than the deeply bound Efimov 
trimers, and they should therefore be easier to study experimentally. Also the lowest trimer energy, 
which most closely follows a free-space Efimov state, could be easier to study in experiment. The trimer 
could first be prepared in the lowest harmonic oscillator level for a/aho <C — 1, and then the scattering 
lenght can be changed via a magnetic field sweep to the region a/aho ~ ~2, where the decay is still 
manageable. 



4 Effective range corrections 

In order to improve the approximation of the channel eigenvalues A„ arising in ( |14[ ), one might want 
to include the influence of the effective range r^. This can be done by making the so-called effective 



range expansion, i.e. using the expansion in Eq. ( 19 1 up to second order. In literature there have been 



various instances in which the effective range is taken into consideration in the determination of the 
three-particle spectrum [331 131] • 

Having an additional interaction parameter at hand, raises the question whether in this regime 
Efimov physics in ultracold gases can still be universal. On the experimental side, there are presently 
also some conflicting observations. In order to understand these issues, we first need to define what we 
mean by the universal Efimov regime. In the above, we have already seen that besides the scattering 
length, we need at least one additional parameter that fixes at least one of the Efimov trimer energies. 
The other ones can then be derived from this position. This information is provided via the three- 
body parameter k*, which can be regarded as a boundary condition for the interaction potentials in 
the regime where the three particles are very close together. It is related to the energy of the lowest 
Efimov trimer via Et = —fi^Kl/m in the infinite scattering length limit, and with this only additional 
parameter the Efimov physics is universal. 

However, some experiments on three-body recombination do not imply universality. While the 
experiments at Rice University did show the universal connection between Efimov states in ''Li at one 
side of the resonance [5], the universality did not extend to the other side of the resonance. The Rice 
group was not alone in this observation: also in Florence [21 on '^^K, in Innsbruck on ^"^'^Cs, and in 
Tokyo [13] on ^Li this non-universality was observed. On the other hand, in three other experiments 
universality across the resonance was observed: in Bar-Ilan it was found that universality does extend 
across a Feshbach resonance in ^Li, and even that the short range physics encapsulated in At* is 
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100 




Fig. 6 Eigenvalues A of depending on R/Ro, as calculated in three different ways. The blue lines represent 
the approximations of the ei geny alues by a truncated Fourier expansion. The black lines and red lines are the 
approximations obtained by (201 and ( |37[ ) respectively. We take for the square well parameters VoRq = —3.2. 
It can be clearly seen that the effective range gives a much better asymptotic behavior. 



independent of the nuclear spin configuration |35l I36j . Furthermore universality was also observed in 
a three-component fermionic ^Li gas in Heidelberg [S] and in State College [TT] . 

Now the question arises whether the effective range paramer rs can explain, in addition to the 
three-body parameter k*, the non-universal behavior of some of the above-mentioned experiments. 
However, it might also be that in a certain regime, the effective range will take over the role of the 
three-body parameter due to its regularising effect of the radial wavefunctions |33l I34| . And recently, it 
has been shown that small experimental deviations from universality 35^ can be explained by effective 
range corrections using an effective field theory [S7] . 

As it has been shown by Fedorov and Jensen [3H] and more recently by Platter et al. [29 , the 



introduction of the effective range in the approximate equation ( 18 1 for the eigenvalues A„ yields the 
relation 



cos 



R , 1 A„ 



2 2i? 



(37) 



To get an idea on the validity of the approximation (37), we have compared the results with the 



eigenvalues obtained by a truncated Fourier expansion for a square well potential with height Vq and 
width i?o, such that VqR^ = —3.2. This expansion is very convenient when one consi ders the operator 
J-^, which depends on the hyperradius R and is defined via the eigenvalue equation (14) as follows 



92 



,{R, a) + 2R^V{V2R sin a) 



i{R,a) 



V3 



|i7r-a| 



4>n{R, ct')da' 



The reason for this convenience is that in the basis of Fourier expansion functions 

2 



: sin(2na), 



(38) 



(39) 



both the differential operator and the integral operator are diagonal. The matrix elements are then 
given by 

(-F«)™„ = An^S^r. + 2R^ (^1 + -^^ sin(2nf )^ V„„(i?), (40) 



where 



Vmn{R) 



Tr/2 



V{V2R sin a) sin(2nQ;) sin(2mQ;)dQ;. 



(41) 



12 




In Figjm we compare the numerically calculated eigenvalues A„(i?) of J-^ with results obtained from 
Eq. ( |20y ^ which depends on the scattering length only, and with the results from solving Eq. (37), 
which depends both on scattering length and effective range. Here the values for a and are taken 
from the effective range expansion of the square well potential. It can be seen that asymptotically, 
that is for large hyperradius, the inclusion of the effective range in the model yields results that are 
considerably better than if one only takes the scattering length into account. For relatively small values 
of the hyperradius, the predictions obtained by the Fourier expansion and this simple model start to 
deviate. Hence, we expect that for short hyperradii, additional details of the interaction potential 
become important, i.e. higher order expansion coefficients beyond the effective range approximation 
have to be taken into account. 



For Ts < and a fixed value for a, the lowest solution Ao to Eq. (37 1 can be defined as a smooth 
function of i? for i? e (0, oo). Moreover, Ao(-R) = 0{R) as i? — 0. Therefore, if we use this approxima- 
tion in the hyperradial equation, we can find regular solutions instead of rapidly oscillating ones. This 
could, at least in principle, remove the need for a renormalization boundary condition in terms of k.^ for 
i? — ^ 0. In the absence of a harmonic oscillator, and in the limits of small energy and r = —r^j a — > 0, 
we will see that the solutions to the hyperradial equation are identical to those when the effective 
range is not taken into account, except that the value of is determined by by boundary layer 
behavior for small values of the hyperradius. Such a relation between and rg has been established 
before [33 El. 



Let us, therefore, study the equation (37 1 more closely. We have plotted the solution Ao for the 
situation a = —20rs in Fig. [7j to give a characteristic example. In the region — ^ i? < a, the value 
of Ao satisfying (37 1 is close to that obtained by solving (20 1. Similarly, when R ^ — Ts, the solution 
Ao to (37) will be close to fj,{—R/rs), where ^{—R/rg) is the solution to 



cos 



V3 



-1/2 



sin(Ay2|)=_y2A-V2 3in(^i/2|) 1, 



'2R' 



(42) 



Around R — r^, there is a transition from a boundary layer behavior similar to /i to behavior obtained 
by ignoring the effective range ( 20 ) . 



To find the low-energy behavior of the system, which is needed to find a proper relationship between 
Ts and [331 134] . we have to turn back to the hyperradial equation (21). Following from the above 
considerations, the solution /o that decays exponentially for large hyperradius is to leading order in 
— r^/a, and for R ^ —r^, equal to /°"*, the exponentially decaying solution to the hyperradial equation 
obtained when the effective range is ignored. This outer solution /°"* will have log-periodic behavior 
as i? 

- V^sin(solog(Ki?) +ao), (43) 

for some k. The solution /o will be similar for —rg ^ i? <C a, and will therefore for low enough energies 
also show log-periodic behavior. However, there is a boundary layer for R <^ —rs, which makes that we 
can demand the solution fo{R) to be regular for i? — )• 0, that is, such that fo{R) ~ VR for i? — >• 0. This 
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effectively determines the value of k. To analyze the boundary layer, we introduce the new variable 
p = —R/ts and the function /o, /o(p) — foi^pfs)- The hyperradial equation for /□ reads 



, Ao(p) - 1/4 ^ rt 



foip)=0. (44) 



and we demand that fo{p) behaves like ^/p for /? — > 0. We consider now the case of a negligible harmonic 
oscillator potential, i.e. Oho — oo, and take the low-energy limit Er^ — > and a diverging scattering 
length. The solution will then be equal to which satisfies 



d' ^ A^(p)-l/4 



= 0. (45) 



For p — >• oo, 

lim (p) ^ sin(so log(p) + <^o) ■ (46) 

p— >-oo 

We determine the phase 4*0 numericahy, by choosing a boundary condition 



fM\^^MM^ (47) 
Vp J p Vp' 



for small enough values of p, and by integrating the equation. The boundary condition ensures that 
we pick up the right solution. The result is 0o = log(1.0918) + ag. In order to find what value of 
(= K in this limit) is induced by r^, we have to match the inner and the outer solutions in the regime 
where — R a. This yields that —K^,rs — 1.0918. 

Remarkably, this numerical value for — K*rs differs from the values obtained in Ref. j33j and |34| . 
which are —K^Vg « 5 and — K^rg « 5.3062 respectively. It is unclear where the difference comes from, 
however, both papers use models different from the hyperspherical formalism that we used here, where 
we base our treatment of the effective range parameter on the descriptions in Refs. ^28^ and [22]. More 
research is needed to resolve this discrepancy. 

Now we are in a good position to compare our model with earlier work on the Efimov effect in a 
harmonic potential. In the work of Th0gersen, Fedorov, and Jensen, the influence of the effective range 
is studied on the universal Efimov effect |25j . They considered a system with a Gaussian interaction 
potential, given by 

V{r)^VoeM~ryRl), (48) 

where Vq and Rq are constants. The range Rq has been chosen to equal Rq = aho/3965. 

Although the main purpose of Ref. [2S] is to show the influence of the effective range, their model 
contains a harmonic trap as well. In order to make a proper comparison with our model, we take 
3.7. The results are shown in Figs. [8] and Fig. |9) and can be directly compared with Fig. 1 
in Ref. . We observe that the numerical results are reproduced quite closely and we can assign the 
numerically found bound states to the different hyperspherical channels. We also observed that if we 
include the above corrections that come from the effective range, the comparison is better than results 
obtained without these corrections. 



5 Conclusion 

We studied the Efimov effect in a harmonic potential by using the hyperspherical formulation. In the 
first part of this paper, we studied the universal regime by only taking a large scattering length into 
account, and by making the adiabatic hyperspherical approximation in that we ignored the coupling 
between different hyperspherical channels. The model thus obtained is simple in the sense that the 
binding energies of the trimers result from an ordinary differential equation, with an appropriate 
boundary condition for the lowest channel encoded by the three-body parameter k*. Considering the 
lowest channel, we formulate an equation for a phase variable, the exponent of which is the radius 
in a spherical plot describing a surface of values in {a~^ ,K,a^^) -space (inverse scattering length, 
wavenumber and inverse oscillator length) for which bound states exists. Other such surfaces are 
obtained by multiplying the distance to the origin by a factor e'^/*'", and there is an accumulation 
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sgn(a)|aho/a|''"* 



Fig. 8 Trimer s olut ions in a harmonic oscillator, for potential parameters corresponding to the paper of 
Th0gersen et al. [25]. Shown are the solutions above the three-particle threshold. The zeroth (first) channel 
solutions are given by black-solid (red-dashed) lines. The graph can be directly compared to Fig. 1 in this 
paper. 
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Fig. 9 Trimer s olut ions in a harmonic oscillator, for potential parameters corresponding to the paper of 
Th0gersen et al. [25]. Shown are the solutions below the three-particle threshold. The red line indicates a 
solution when the effective range is taken into account. Also indicated is the dimer-atom threshold, which 
includes effective range corrections. The graph can be directly compared to Fig. 1 in this paper. 



point of surfaces near the origin, analogous to the Efimov effect in free space. In the limit aho — >■ oo, 
the Efimov effect in free space is recovered. 

Interesting differences with the situation of three particles in free space can be observed when we 
focus at a fixed value of aho- For parameter combinations where in free-space bound states would exist, 
the situation is very similar, except that the accumulation of energies only happens up to a typical 
energy of hiu. The rest of the Efimov states are "pushed out" into the trimer continuum (which is of 
course only a continuum in free space). Higher hyperspherical channels give rise to additional bound 
states, and by taking those into account as well we obtain results very close to numerical simulations 
of three-particle systems |25| . Better qualitative results are obtained by approximating the coupling 
between the hyperspherical channels, which lifts the degeneracies and creates avoided crossings. The 
lifetime of the trimers is highest for states with a continuum character, i.e. in the region which we 
would denote as trimer continuum in free space. The reason for this is the small overlap of these states 
with deeply bound states. Experimental study of Efimov-like states, however, could be feasible close to 
the transition from this continuum region to the trimer region, where the linewidth is still quite small. 
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For a more general Efimov model that allows for deviations from the universal regime, a natural step 
is to improve on the approximation of the hyperangular problem by taking into account the effective 
range. By doing this, we were able to find a relation between and the effective range rg < if 
a ^ — Ts 3> -Roi where i?o is a typical range of the potential. We obtain the same qualitative behavior 
as Ref. |33| and [33] where different descriptions were used, however, the actual value of kTs is also 
very different. It appears that the actual value for this relationship depends strongly on the choice 
of the model and the corresponding approximations to describe the three-particle physics, and more 
research is needed to understand this discrepancy. In a different parameter regime for a, and Rq, 
such a relationship between and is not necessary. Then we have in addition to , the three-body 
parameter in the boundary condition for the hyperradial problem to our disposal. By comparing to 
numerical calculations we find that this description indeed improves the predicted trimer energies. 
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